Bulk Tunneling at Integer Quantum Hall Transitions 
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The tunneling into the bulk of a 2D electron system (2DES) in strong magnetic field is studied near 
the integer quantum Hall transitions. We present a nonperturbative calculation of the tunneling den- 
sity of states (TDOS) for both Coulomb and short-ranged electron-electron interactions. In the case 
of Coulomb interaction, the TDOS exhibits a 2D quantum Coulomb gap behavior, v(e) = CQ|e|/e 4 , 
with Cq a nonuniversal coefficient of quantum mechanical origin. For short-ranged interactions, we 
find that the TDOS at low bias follows u(e)/u(0) = 1 + (|e|/eo) 7 , where 7 is a universal exponent 
determined by the scaling dimension of short-ranged interactions. 



PACS numbers: 73.50.Jt, 05.30.-d, 74.20.-z 

The integer quantum Hall transition (IQHT) refers 
to the continuous, zero-temperature phase transition be- 
tween two consecutive integer quantum Hall states in a 
2DES Jl],^| . It happens when the Fermi energy of the 
2DES moves across a critical energy located near the 
center of each disorder-broadened Landau level. On both 
sides of the transition, while the Hall conductivity a xy is 
integer-quantized, the dissipative conductivity a xx van- 
ishes at low temperatures. At the transition, however, 
both a xx and a xy remain finite, supporting a critical, 
conducting state in two dimensions. 

Theoretically, the IQHT has been studied extensively 
in terms of its important, noninteracting analog in which 
electron-electron interactions are ignored §]. The ex- 
tend to which this disordered, free electron model de- 
scribes the IQHT in real materials depends on the effects 
of electronic interactions. Recently, the stability of the 
noninteracting theory has been analyzed by calculations 
of the scaling dimensions of the interactions || . It was 
found that short-ranged interactions are irrelevant in the 
renormalization group sense, but the long-ranged 1/r- 
Coulomb interaction is a relevant perturbation. In the 
latter case, the noninteracting theory becomes unstable, 
and the universality class (e.g. the critical exponents) of 
the transition is expected to change based on the stan- 
dard theory of critical phenomena. The experimentally 
measured value of the dynamical scaling exponent, z = 1, 
is consistent with Coulomb interaction being relevant at 
the transitions [Q. However, the precise role played by 
Coulomb interaction and the mechanism by which it gov- 
erns the scaling behavior have not been understood. 

The tunneling density of states (TDOS) is a simple and 
useful probe of the nature and the effects of electronic in- 
teractions. Recent tunneling experiments in the integer 
quantum Hall regime discovered, remarkably, that the 
TDOS vanishes linearly on approaching the Fermi en- 
ergy @- Since at the transition, a xx is finite and the 
localization length £ is very large, the linear suppression 
of the TDOS at criticality is expected to have a differ- 
ent origin than the 2D classical Coulomb gap behavior 
deep in the insulating phases. The results of numerical 
calculations of the TDOS using Hartree-Fock approxi- 



mation (HFA) of Coulomb interaction also show a linear 
Coulomb gap at all Fermi energies in the lowest Landau 
level ||[7| . Since HFA does not include screening of the 
exchange interaction, the correct behavior of the TDOS 
remained unclear, especially at the transition. 

In this paper, we present a nonperturbative calculation 
of the TDOS at the IQHT, taking into account the dy- 
namical screening of Coulomb interaction by the diffusive 
electrons. We show that the TDOS is given by, 



v(e) = C Q \e\/e 



(1) 



at low energy e (bias). We shall refer to Eq. (|l|) as the 
2D quantum Coulomb gap behavior. The coefficient Cq 
is not a universal number as in the 2D classical Coulomb 
gap expression || , but rather a quantity of quantum me- 
chanical origin. In general, it depends on magnetic field 
and the microscopic details of the sample such as the 
mobility. We find, at the IQHT, 
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where g = 2ir 2 a c with a c the critical conductivity in units 
of e 2 /h and A = (A;/Ob) 2 /(27t 2 with as the Bohr 
radius and / the zero-field mean free path. $ is the error 
function. For large e, v(e) crosses over to the perturba- 
tive diagrammatic result in strong magnetic fields ||[l0) . 
We also study the case of short-ranged interactions corre- 
sponding to experimental situations when the presence of 
nearby ground planes or the tunneling electrodes them- 
selves screen the long-ranged Coulomb interaction be- 
tween electrons are large distances. In this case, we find 
that the zero-bias TDOS is suppressed from the nonin- 
teracting value but remains finite. The corrections at fi- 
nite e follow a universal power law: \v{e) — v(0)]/v(fi) = 
(|e|/eo) 7 with 7 = x/z, where x ~ 0.65 is the scaling 
dimension of short-ranged interactions and z = 2 [jl . 
We begin with the partition function, 
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where the imaginary time action S — Sq + Si, 
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Here, if) and i\> are the electron Grassmann fields, <p is 
a scalar potential coupled to electron density, and A ai 
a = x, y is the external vector potential, e a pd a Ap = Bz. 
Integrating out cj> in Eq. (||) gives rise to the electron- 
electron interaction U(r — r') in Coulomb gauge. Quench 
average over the short-ranged impurity potential V is 
represented in Eq. (^) by integrating over a Gaussian 
distribution P[V] using the standard replica method. 

The TDOS can be obtained from the impurity aver- 
aged Green's function G(r) = (ip(r,T)tp(r,0)), 
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In the imaginary time path integral, G(r) is given by 
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Consider a U(l) gauge rotation, ip — > tpe 10 , %j) — > e~ lB 
with 9 = 6(r,r). In the rotated frame, G(r) has the 
same form as in Eq. (ra), but with a transformed action 



S T = S -i d rdt(jla - j M a M ), 



(7) 



where a M = d^O, jj, 



T,x,y, is the longitudinal U(l) 
gauge field coupled to the fermion 3-current , and 



(8) 



is the source density (current) disturbance due to the tun- 
neling electrons. Because of the latter, the ground state 
develops a nonvanishing charge and current density, and 
the saddle point of the tunneling action S T is shifted from 
that of S in Eq. (^). Minimizing the action, dS/d9 = 0, 
one finds that the induced charge and current density, 
P = (jo) an d J = (j), satisfy the continuity equation, 
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where T(r,t) — d T jo — S(r)[S(t) — 5(t — r)] corresponds 
to the process of injecting an electron at r = and time 
t = and removing it at t — r. Since the critical conduc- 
tivity is finite at the IQHT, the charge spreading is ex- 
pected to be described by (anomalous) diffusion, we have 
J = — D(V — x V)p where D is the field-dependent 
diffusion coefficient and jh — o~ xy /cr xx is the Hall an- 
gle. Inserting this into Eq. (||), we obtain the diffusion 
equation, 



-p-DV 2 p = T{r,t). 



(10) 



Notice that 7^ does not enter because the transverse 
force does not affect the charge spreading in the bulk of 
the sample. Solving Eq. (HO) leads to 
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where T(r,u> n ) = 1 — exp(— zw„t) and to n is the boson 
Matsubara frequency. 

Next we go back to the tunneling action in Eq. (Q) and 
expand the currents around the expectation values: jo = 
P + $jo an( i j = J + <5j . The fluctuating parts now satisfy 
the homogeneous continuity equation: dtSjo + V • Sj = 0. 
For convenience, we choose the unitary gauge by setting 
4>{r, t) — a (r, t). Using Eq. (Q), the tunneling action H] 
for charge spreading becomes, 



S T = S u — i / dtd rpao, (12) 
Jo 

where the unitary gauge action S u is given by 

S u = So + i / dtd rSj a a a (13) 
Jo 
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dtd 2 rd 2 r'a (r)U- 1 a (r'). 
The one-particle Green's function in Eq. (||) becomes, 



G(r) = Z~ 



V[i>, ^}V[9]V[V)P[V^(T)tP(0)t 



(14) 



We now quench average over the impurity potential and 
integrating out the fcrmions for a fixed gauge configura- 
tion @,[l3|. Eq. (ph becomes, 
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G(r) = / V[6]G e (T)e 
where Gg(r) = (ip(T)ip(0))e is given by, 
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with the gauge field effective action, 

e -S ef f[o M ] = / v[iP,$]V[V]P[V}e~ s ^^' a » ] . 



(15) 
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Thus far, all evaluations have been formal and exact, 
and explicitly displayed the structure of the theory in 
the U(l) sector. 

To calculate G(r) in Eq. ( |l5| ) , one has to make system- 
atic approximations to obtain S c g [a M ] and GgM. Follow- 
ing essentially the formalism of Finkelshtein [0-15[ , we 
obtain S e g to quadratic order in a, 



2 



Seff = T 



d 2 qa (q,Lo n )U(q,Lu n )ao(-q,-Lo n ). (18) G(r) = G (t) exp 



Here II is the polarization function, 
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which coincides with the dynamical screened Coulomb 
interaction ||[l0]]- In Eq. ([l9|), dn/djj, is the compress- 
ibility. The bare diffusion constant in the self-consistent 
Born approximation (SCBA) is D = ^r 2 ^ 1 , where 
r c = (2N + l) 1 / 2 /^ with Ib the magnetic length and 
N the Landau-level index. The effects of anomalous dif- 
fusion will be discussed later. Notice that in strong field, 
D is proportional to the field dependent scattering rate 
l/r (£)~ V^Ao(0)._ 

Next we turn to Ge in Eq. (|l6|). In the unitary gauge, 
the important interference effects between the phase of 
the electron wave functions have been accounted for in 
the tunneling action in Eq. (p"5|). The amplitude fluctu- 
ations are small for the slowly varying gauge potentials 
that dominate the dynamically screened Coulomb poten- 
tial in Eq. rtljJ) p3| . This is a unique feature of the slow 
diffusive dynamics of the electrons. We therefore neglect 
the ^-dependence and write Gg(r) « Gq(t). The func- 
tional integral over the gauge potential in Eq. (JlJ) can 
thus be carried out explicitly using the effective action in 
Eq. (pi), leading to 



G(r) = G Q {T)e w[T \ 



(20) 



For Coulomb interaction, U(q) — 2ne 2 /q in Eq. (19) 
Making use of Eq. (Elf), W(r) has the form, 
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where k = 2ire 2 dn/dfi is the inverse screening length at 
the transition. Notice that W(r) has a similar struc- 
ture as the leading correction to the TDOS in the di- 
agrammatic perturbation theory in both zero jl6| and 
strong magnetic field |J. The main contribution to the 
q- integral in Eq. ( pl| ) comes from the region \ui n \/Dn < 
q < {\u n \/Dyi 2 . In this region, the diffusion coefficient 
D is a constant. The anomalous diffusion |H only ap- 
pears in the opposite limit Dq 2 \u> n \. As a result, 
W(t) becomes, in the T = limit, 
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where t s — 1/Dn 2 , and the conductivity a xx = Ddn/d^, 
in units of e 2 /U following Einstein's relation. The up- 
per limit of the integral in Eq. ( p2|) is Tl/tq since one 
can show t s <C tq near Landau level centers. The 
long-time behavior of W is therefore given by W(t) — 
(-l/87r 2 ^ K )hi(r/ro)ln(r/(T s 2 /To)) @. Substituting 
this result into Eq. (pQ), we obtain, 



(23) 



The asymptotic behavior of Go (r) for large r is Go (r) ~ 
—vq/t, where is the corresponding DOS |l3| . In 
the high-field limit of the SCBA, it is well known that 
v Q = (l/2Trl 2 B )(2T )/h and a xx = v Q D ~ (2A^ + 1)/2tt 2 
in the center of the N-th Landau level. After analytical 
continuation to real time [p^[ , we obtain the final result 
for the TDOS defined in Eq. (|), 
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We now discuss the behavior of v(e) in different 
regimes, after making a few remarks, (i) The quantity 
in the exponential in the above equation did not follow 
from an expansion in l/a XXl but rather resulted from the 
leading contribution dominated by the anomalously di- 
vergent ln 2 -term at long-times. Next order corrections 
to the latter are of the order {l/a xx , 1/<j xx } lnt/r . (ii) 
While Eq. (^4|) correctly captures the double-log contri- 
butions, we have assumed that the conductivity a xx does 
not depend on frequency in Eq. (|2^). However, o ' xx can 
be renormalized by localization effects (of leading order 
(l/cza;) \vllotq in the unitary ensemble) and interaction 
effects (of leading order lnwro in strong magnetic field 
|)].|lo|). Thus o xx assumes the frequency- independent 
SCBA value only if — ln(ero) o xx . In this regime, 
and for — ln(ero) ^> ^/a xx , the integral in Eq. (|24| ) gives, 



v(e) = Vq exp 
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ln(| £ |ro)ln(| £ |r 2 /r ) 



(25) 



Expanding the exponential to leading order in l/a xx , 
Eq. ( p5| ) reproduces the high-field diagrammatic pertur- 
bative result of Girvin, et. al. ||. In the case of B = 0, 
such nonperturbative resummation of the perturbative 
double-log divergences was pointed out by Finkclshtcin 
[ 14 1 , and recently reexamined using different approaches 
|llf,RJJl9). The result of Eq. @ can be regarded as an 
extension of the latter to high magnetic fields. 

(iii) As stated in (ii), at low-bias, i.e. for — ln(er ) ^> 
a xx , the localization and interaction effects lead to, in 
general, a frequency dependent conductivity <t xx (uj) in 
Eq. (IH). However, at the IQHT, the critical conduc- 
tivity a c is finite and of the order of e 2 /h p| j7|, p0| , pl| ■ 
Thus, the validity of our analysis, i.e. the structure of 
the double-log divergence at long-times, can be extended 
into the regime of small e, provided that a xx in Eq. ( |24| ) is 
replaced by the critical conductivity a c . It is important 
to note that because of the double-log term, the expo- 
nential factor in Eq. ( |24| ) converges very fast such that 
the TDOS becomes analytic in small |e|. Expanding to 
leading order in |e|, 
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Performing this integral, and use the fact that the com- 
pressibility is only weakly renormalized, i.e. dn/dn ~ v , 
we obtain the results given in Eqs. (|l|) and (g) with 
A = t s /tq. We thus conclude that for the 1/r-Coulomb 
interaction, the bulk TDOS at IQHT exhibits a linear 
Coulomb gap behavior of quantum mechanical origin, i.e. 
the 2D quantum Coulomb gap. As a consequence, the 
level spacing at the transition becomes AE ~ where 
L is the size of the system, leading to a dynamical scaling 
exponent z = 1. This behavior of the TDOS is qualita- 
tively different from those obtained in the clean case ]22| ] 
and in a weak magnetic field . 

We emphasize that the linear quantum Coulomb gap 
behavior results from the combined effects of (i) two- 
dimensionality, (ii) 1/r — Coulomb potential, and (iii) 
quantum diffusion, i.e. a finite conductivity at T — 0. ft 
pertains therefore to other metal-insulator transitions in 
2D amorphous electron systems, provided that the crit- 
ical conductivity is finite An example is the re- 
cently discovered 2D B = metal-insulator transition 
p4| , although the asymptotic low temperature behavior 
of the critical conductivity extracted from the experimen- 
tal data in this case is still controversial. 

Finally, we turn to short-ranged interactions, ft has 
been shown interacting potentials f jr v that decay faster 
than 1/r 2 are irrelevant and scale to zero with scaling 
dimensions x = p— 2 for 2 < p < 2 + X4 S , x^s ~ 0.65, and 
x = X4 S for p > 2 + X4 S at the IQHT jj) . For simplicity, 
we will consider U(r — r') = uS(r — r'), thus U(q) = u in 
the screened interaction in Eq. (]l9|) . Inserting the latter 
into Eq. (21), the q integral no longer leads to the log- 
divergent term in to as in Eq. (|22| ) due to the short range 
nature of the interaction. Instead, we find |l5| ], 



W„(t) = - 
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where a is a nonuniversal quantity dependent on the in- 
teracting strength. It is given by 
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with A = up 



C = 1/2 + 1/(2 - 3r?/2), and n the 
anomalous diffusion exponent p8[ . Once again in this 
case, the critical conductivity a c is finite at the transi- 
tion such that Eqs. ( pi7j ) and ( p8| ) represent the leading 
contribution to W sy {t) in the asymptotic limit. If the 
interaction u were a marginal perturbation, then one 
would obtain, as in the Luttinger liquid case, W sr = 
-aln(r/r ), G sr (r) - T ~^+ a \ and v{e) ~ v \e\ a . How- 
ever, this is not true, because u is an irrelevant per- 
turbation and the effective interaction scales to zero ac- 
cording to itoff ~ wjj x l z with z = 2 at the noninteract- 
ing fixed point pJ. As a result, a obeys a scaling form 
a(u,ui) = i(uo) 1 ' 2 ). The fact that A(u -> 0,w) = im- 
plies, together with Eq. (p8|), the leading scaling behav- 
ior a ~ AX(ujt ) x / z with 4 = C/47r 2 cr c . Substituting this 



into Eq. @, one finds that W sr (r) = ^^[(to/t^-I] 
where j = x/ z ~ 0.32. That W sr (r) converges in the 
large-r limit is a consequence of the interaction being 
irrelevant, i.e. 7 > 0. Thus, we find for large r, 



G sr (r) = ^(0)-exp 

T 



AX 

1 



(29) 



where u(0) = vq exp(— AA/7) < i/q. After analytic con- 
tinuation, the TDOS at small e is given by, 



v„(e) = i/(0) 



1 + 











(30) 



where e = t _1 (AA/7) _1//7 . This result leads to sev- 
eral interesting predictions: (a) For short-ranged inter- 
actions, the TDOS is finite and nonuniversal at zero bias 
i/ SI (s = 0) = v(0) ^ 0. (b) Since v(0) < v for large 
A, interactions still lead to strong density of states sup- 
pression at low bias, (c) The plus sign in Eq. ( ^o|) indi- 
cates that the TDOS increases with bias e according to 
a universal power law with an initial cusp singularity for 
our value of 7. These predictions can, in principle, be 
tested experimentally by deliberately screening out the 
long-ranged Coulomb interaction using metallic gates. 

The authors thank R. Ashoori, M. P. A. Fisher, S. M. 
Girvin, I. A. Gruzberg, A. Kamenev, Y. B. Kim, D.-H. 
Lee, and N. Read for useful discussions. This work was 
supported in part by NSF Grant No. PHY94-07194 (at 
ITP), DOE Grant No. DE-FG02-99ER45747, and an 
award from Research Corporation. 



[1] For reviews see, The Quantum Hall Effect, edited by R.E. 
Prange and S. M. Girvin (Springer- Verlag, New York, 
1990); Perspectives in the Quantum Hall Effects, edited 
by S. Das Sarma and A. Pinczuk (John Wiley & Sons, 
New York, 1997). 

[2] B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 

[3] D-H Lee and Z. Wang, Phys. Rev. Lett. 76, 4014 (1996). 

[4] L. W. Engel, et. al. Phys. Rev. Lett. 71, 2368 (1993). 

[5] H. B. Chan, et. al. Phys. Rev. Lett. 79, 2867 (1997). 

[6] S.-R. E. Yang and A. H. MacDonald, Phys. Rev. Lett. 
70, 4110 (1993); S.-R. E. Yang, A. H. MacDoanld and 
B. Huckestein, Phys. Rev. Lett. 74, 3229 (1995). 

[7] B. Huckestein and M. Backhaus, Phys. Rev. Lett. 82, 
5100 (1999). 

[8] A. L. Efros and B. I. Shklovskii, J. Phys. C8, L49 (1975). 
[9] S. M. Girvin, M. Jonson, and P. A. Lee, Phys. Rev. B26, 
1651 (1982). 

[10] A. Honghton, J. R. Senna, and S. C. Ying, Phys. Rev. 

B25, 6468 (1982). 
[11] L. S. Levitov and A. V. Shytov, Sov. Phys. JETP Lett. 

66, 214 (1997). 

[12] Y. B. Kim and X. G. Wen, Phys. Rev. B50, 8078 (1994). 

[13] S. Xiong and Z. Wang, to be published. 

[14] A. M. Finkelshtein, Sov. Phys. JETP 57, 97 (1983). 



4 



[15] A. Kamenev and A. Andreev, cond-mat/9810191 



[16] B. L. Altshuler, A. G. Aronov, and P. A. Lee, Phys. Rev. 
Lett. 44, 1288 (1980). 

[17] For a general bare potential U(r) ~ l/r p , one can show 
jl3j that for p < 2, diffusive screening dominates and the 
long-time behavior of W(t) has the double-log form. The 
case of p > 2 is similar to the short-ranged interaction 
case discussed below. 

[18] J. T. Chalker and G. J. Daniell, Phys. Rev. Lett. 61, 
593 (1998). 

[19] P. Kopietz, Phys. Rev. Lett. 81, 2120 (1998). 
[20] Huo, R. E. Hetzel, and R. N. Bhatt, Phys. Rev. Lett. 
70, 481 (1993) 

[21] Z. Wang, B. Jovanovic, D-H Lee, Phys. Rev. Lett. 77, 
4426 (1996). 

[22] S. He, P. M. Platzman, and B. I. Halperin, Phys. Rev. 

Lett. 71, 777 (1993). 
[23] I. L. Aleiner, H. U. Baranger, and L. I. Glazman, Phys. 

Rev. Lett. 74, 3435 (1995). 
[24] S. V. Kravchenko, et. al. Phys. Rev. B50, 8039 (1994); 

ibid, 51, 7038 (1995). 



5 



